STM induced hydrogen desorption via a hole resonance 
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We report STM-induced desorption of H from Si(100)-H(2x 1) at negative sample bias. The 
desorption rate exhibits a power-law dependence on current and a maximum desorption rate at 
— 7 V. The desorption is explained by vibrational heating of H due to inelastic scattering of tunneling 
holes with the Si-H 5a hole resonance. The dependence of desorption rate on current and bias is 
analyzed using a novel approach for calculating inelastic scattering, which includes the effect of the 
electric field between tip and sample. We show that the maximum desorption rate at —7 V is due 
to a maximum fraction of inelastically scattered electrons at the onset of the field emission regime. 
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Single atom manipulation with a scanning tunneling 
microscope (STM) has been reported for several systems 
and a variety of physical mechanisms has been proposed 
to account for such manipulation |l],H . Among the works 
relevant for the present Letter we mention single hydro- 
gen atom desorption from the Si(100)-H(2xl) surface [3|, 
and dissociation of single O2 molecules on Pt(lll) [4|. 
These manipulations were performed at positive sam- 
ple bias, and the underlying microscopic mechanism has 
been related to vibrational heating by inelastic scattering 
of tunneling electrons with an electron resonance on the 
sample. There have been theoretical predictions that a 
related mechanism may operate at negative sample bias 
[gj , which involve inelastic scattering of a tunneling hole 
with a hole resonance on the sample. However, to our 
knowledge there has been no experimental confirmation 
of such a mechanism, probably because high tunnel cur- 
rents and sample biases are needed to obtain high inelas- 
tic scattering rates with low-lying hole resonances (see 
Fig. 1). 

In this Letter, we present evidence for a desorption 
mechanism involving a hole resonance, for STM-induced 
hydrogen desorption from the monohydride Si(100)- 
H(2xl) surface in ultra-high vacuum Whereas in 
previous studies J3|,|],0, hydrogen desorption has been 
studied at positive sample bias, here we report hydrogen 
desorption at negative bias. The desorption process is 
modelled by vibrational heating of hydrogen caused by 
inelastic scattering of tunneling holes with the Si-H 5a 
hole resonance. The inelastic scattering rates are calcu- 
lated using a novel method based on first principles elec- 
tronic structure theory, and desorption rates obtained by 
solving the Pauli master equation for a truncated har- 
monic potential well j| . With this model we obtain des- 
orption rates as function of tunnel current and sample 
bias that are in agreement with the experimental data. 
We find a maximum desorption rate at sample bias — 7 V, 
which coincides with the onset of the field emission or 
Fowler- Nordheim regime ||[l0| . We show that this results 
from a competition between polarization of the hole res- 



onance, which increases the fraction of inelastic scattered 
electrons, and domination of Fermi-level contributions to 
the tunnel-current in the field emission regime. 

The experiments were performed on n-type Si(100) 
(N D = 1 x 10 18 cm" 3 ) samples using a JEOL JSTM- 
4000XV microscope at a base pressure of 1 x 10 -10 torr. 
Atomic hydrogen was absorbed on the clean Si(100)- 
(2x1) reconstructed surface to obtain the monohydride 
(2x1) phase, in a manner identical to previous reports 
0, and W tips were used. The desorption experiments 
were carried out by scanning the STM tip at speed, s, 
sample bias, Vb, and tunnel current, /, and subsequently 
imaging the affected area to determine the number of 
Si sites where desorption occurred. Figure 2 shows a 
typical example of desorption at negative sample bias of 
hydrogen along a single dimer row. We first measured 
the dependence of the desorption rate, R, on the tun- 
nel current for sample biases of Vb — —7 V and —5 V. 
The results are shown in Fig. 3(inset). For both biases 
there is a power-law dependence of the desorption rate 
upon current, R = Rq(I / Idcs) a , with exponent a ~ 6. 
In this equation I^cs is the tunnel current that gives 
rise to a fixed desorption rate Rq. In order to find the 
voltage dependence of the desorption rate, we measure 
Idee{Vh) with Vb in the range —10 V to —4 V, as shown 
in Fig. 3. The measurements were obtained by scanning 
30 nm along a single dimer row at s = 2 nm/s, and I^es is 
defined as the current that gives rise to desorption of 50% 
of the hydrogen along such a scan line, which corresponds 
to R = 4 s -1 . 

The main feature of Fig. 3 is that significantly higher 
bias and tunnel currents are required at negative bias, 
compared to the positive bias case ||. Whereas ides(Vb) 
decreases monotonically at positive Vb |§, it displays a 
minimum at negative Vb, increasing for Vb < —7 V. Ex- 
perimentally, the lithographic resolution at negative bias 
is comparable to that at positive bias. At the highest cur- 
rent levels there is a tendency for the tip resolution to be 
affected by the lithography process. This limits the volt- 
age range over which Ides(Vb) can be readily measured, 
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and may explain why no detailed investigation of desorp- 
tion at negative sample bias has been made previously. 
We note that the behaviour of /des(Vb) at negative sam- 
ple bias has been confirmed independently on a separate 
STM system 0. 

At negative sample bias, electrons accelerated across 
the tunnel gap impinge on the tip, so the possibility of 
direct excitation by collision can be ruled out as a des- 
orption mechanism. In contrast, at positive sample bias, 
direct excitation of the Si-H bond is believed to play the 
dominant role for Vb > 4 V H|). Shen et al. ||| have 
proposed vibrational heating of hydrogen [fl2"[ by tunnel- 
ing electrons scattering inelastically with the Si-H 6cr* 
electron resonance as a desorption mechanism for pos- 
itive sample bias and Vt> < 4 V. A similar mechanism 
can function at negative sample bias, since an electron 
tunneling from the sample to the tip may excite the Si-H 
5<T hole resonance, which upon deexcitation can trans- 
fer energy to the hydrogen atom This process can 
be viewed as inelastic scattering of a tunneling hole with 
the Si-H 5(J hole resonance, as illustrated in Fig. 1. A 
characteristic feature of vibrational heating by inelastic 
scattering is a power-law dependence of the desorption 
rate on current |l3| ] in agreement with the experimental 
observations. 

To make a quantitative analysis of the dependence of 
the desorption rate on tunnel current and sample bias, 
we develop a method for obtaining inelastic scattering 
rates from first principles electronic structure calcula- 
tions. Since the inelastic scattering rate depends on the 
tip-sample distance, d, while the measured quantities are 
I and Vb, we first calculate the relation between /, Vb, 
and d. For this purpose we use a high voltage exten- 
sion of the work by Tersoff and Hamann pi} , that in- 
cludes the effect of the electric field between tip and 
sample 15 . The tunnel current is obtained from the 



local density of states(DOS) of the Si(100)-H(2xl) sur- 
face, p(r, e, E) = |Vv(r, E)\ 2 S(e - e^), at tip position 
r, and wave functions ip^ for electrons with energy 
are calculated in an external field E. The electric field 
is modelled by a planar electric field outside the surface 
and its magnitude determined from Vb and d. The tunnel 
current is given by |l5fl 



= 2fi wK (e) 



\p(d + R w ,e,E)de, 



(1) 



where distances are in bohrs, current in Amperes and all 
electron energies are relative to the sample Fermi energy, 
£f- In this equation k(e) = ^J2m((f)t + eVb — £)/fo is the 
wave function inverse decay length, <j) t — 4.5 eV the work 
function of the W tip |lq] , i? w = 3 a.u. the atomic ra- 
dius of W, and the normalization constant C w = 0.007i? 2 
Amperes x bohr is obtained from a calculation of a model 
W tip (L§. 

We test the method by calculating the STM corru- 
gation across a single missing H defect. The electronic 
structure calculations are based on density-functional 



theory [ |17t|lq ] within the generalized-gradient approxi- 
mation(GGA) JljJ using 20 Ry plane-wave basis sets. 
The Si(100)-H(2x 1) surface is represented by a 12 layer 
(2x1) slab, and we use a c(4 x 4) slab to calculate the cor- 
rugation of a single missing H defect. Ultra-soft pseudo- 
potentials (2^] are used for both H and Si. Figure 2(b) 
shows the result compared with the measured corruga- 
tion perpendicular to a desorption line. The theoretical 
corrugation compares well with the experimental data 
when the electric field between tip and sample is included 
(dashed line) but not otherwise (dotted line). 

We next calculate the tunnel current of inelastically 
scattered electrons. Electrons may scatter inelastically 
due to long-range electrostatic interactions with a vibra- 
tional transition dipolc moment, or due to local inelastic 
scattering with an electronic resonance |2l]]. Only a lo- 
cal interaction can cause atomic scale desorption, so in 
the following we restrict the analysis to resonance cou- 
pling. The resonance scattering event can be described 
by a local polaron model with a linear electron-phonon 
coupling A, and we use a harmonic approximation for 
the Si-H bond potential with frequency ujq. To calculate 
the inelastic current, /„, which causes transitions from 
vibrational level to n, we combine previous models of 
inelastic-scattering [j5|]22]] with the Tersoff-Hamann model 
of STM tunneling between a surface and a W tip Jl4],[l5| . 
In the limit of a broad resonance, A ^ A,cjo, we obtain 
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p n (d + R w ,e,E)de, (2) 



where distances are in a.u. and /„ in Amperes. We de- 
fine a dimensionless parameter, K = irX 2 p s /A, where 
Ps = J2fj, d{ £ 5<? — e n) is the average DOS around the res- 
onance energy £5^, and A is the resonance width. The 
weighted local DOS, 

is weighted with, (K\ (5a\p) | 2 )™, where |5cr) and \p) are 
the resonance state and sample eigenstates, respectively, 
and ffj, = I (5er|^) \ 2 / (x + | (5cr\p) 2 ) is the fraction of elec- 
trons which tunnel from the tip to state p, via the 5a 
resonance. The parameter x determines the fraction of 
electrons which does not tunnel via the resonance state, 
and we have estimated x « 0.1 {— 0.25 x max^KScrl/i)! 2 ). 
We note that the calculations are quite insensitive to the 
value of x, since ~ 1 for e M ~ and / M ~ \(5cr\p)\ 2 / x 
for le^ — £53. 1 3> A, thus x merely damps contributions 
from eigenstates not in resonance. 

We next calculate the parameters entering Eq. (2) for 
inelastic scattering with the Si-H 5cr resonance. From a 
frozen phonon calculation we obtain Hujq = 0.26 eV. Us- 
ing the ground state energy of a free H atom we calculate 
Edcs — 3-36 eV. To find (5cr|/i) we project the electronic 
eigenstates of the slab calculation onto the 5tr wave func- 
tion of a Si-H molecule, and the solid line in Fig. 4 shows 
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the partial DOS n 5(T (e) = \(5a\n) \ 2 8{e - e„). We 
find that the 5cr resonance is centered at £5^ = —5.3 eV 
relative to the Fermi level of an n-type sample with 
an average width of A = 0.6 eV. Crosses in Fig. 4 
show p s (e M ) w n(e M )/|(5(r|/i}| 2 from which we estimate 
p s k 1.2 eV -1 . The electron-phonon coupling is given by 
A = e' 5a y/h/2Mu!o, where M is the hydrogen mass and 
e' 5a . the derivative of £50- with respect to the Si-H bond 
length z H . Our calculations show that £50- varies almost 
linearly with z in the range 1.25 A < z < 2.25 A with a 
slope e' 5c = 2.3 ± 0.1 eV/A, and hence A re 0.20 eV. 

From the inelastic currents we calculate desorption 
rates by solving the Pauli master equation for the tran- 
sitions among the various levels of the oscillator Q. We 
include contributions from I\, I2, and I3 24 1 , and vibra- 
tional relaxations are described by a current independent 
relaxation rate, 7 = 1 x 10 8 s _1 [||J2^]. We assume that 
desorption occurs when the energy of the H atom exceeds 
the desorption energy E^ cs = 3.36 eV, corresponding to 
a truncated harmonic potential well with 13 levels. The 
solid lines in Fig. 3(inset) show the calculated desorp- 
tion rate as function of tunnel current for Vb = —5 V 
and —7 V. The agreement with the experimental data is 
excellent. The desorption rate follows a powerlaw with 
exponent a re 6.5. An analysis of the theoretical calcula- 
tion shows that the contribution from I2 dominates the 
desorption rate, and therefore a re N/2, where N = 13 
is the number of levels in the truncated harmonic po- 
tential well. This contrasts with the low positive bias 
case, Vb < 4 V, where I\ dominates and a re 13 [||. 
This difference is due to the different lifetimes of the 5a 
and 6cr* resonances. Since the 6a* resonance is more 
shortlived(A = leV), the energy transfer in each inelas- 
tic scattering event is smaler and I± dominates in this 
case j23|. Furthermore, electrons with energy < 4 eV 
only sample the low energy tail of the Si-H 6a* reso- 
nance, and p\ is relatively broader than p n for n > 1, see 
Fig. 1. 

We next calculate Jdes(Vb) and the results are shown 
in Fig. 3, and we note that the minimum at —7 V is ac- 
curately reproduced by the theoretical calculation. This 
minimum does not coincide with the resonance energy 
£5<j, but with the onset of the field emission regime, as 
illustrated by I — V characteristics for this system at a 
constant field of 0.8 V/A (dashed line). In field 0.8 V/A 
surface states are bandbent by ~ leV |L5| and therefore 
£ 5cr re 6 eV. Thus, in the range [— 6.5V,0V] the bias de- 
pendence of Ides is mainly determined by the shape of 
the resonance wavefunction. Below —6.5V the bias de- 
pendence is a result of competition between two different 
effects: polarization of the hole resonance and increas- 
ing Fermi-level conduction. The electric field polarizes 
the surface, displacing wave functions towards the tip. 
Since the 5a orbital is more polarized than bulk states, 
its overlap with the tip state increases relative to other 
states, and the fraction of electrons causing double vi- 
brational excitations, fi = I2/I, increases with electric 
field, as observed in the voltage range [— 7V,— 6.5V]. On 



the other hand, the growing tip-sample distance with in- 
creasing magnitude of the negative bias will increase the 
tunnel barrier of the 5a state relative to states near the 
Fermi level, thus decreases with tip-sample distance. 
This effect becomes dominant at the onset of the field 
emission regime, where the conduction of states near the 
Fermi level is independent of the tip-sample distance, and 
f% starts to decrease when Vb < — 7 V. 

In conclusion, we have presented experimental mea- 
surements of the voltage and current dependent variation 
of the hydrogen desorption rate from the Si(100)-H(2x 1) 
surface at negative bias conditions. Based on a novel first 
principles theory of inelastic scattering, we have shown 
that the desorption is caused by vibrational heating of 
the H atom due to inelastic-scattering with the Si-H 5a 
hole resonance. 
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FIG. 1. (a) Inelastic tunneling of a hole into an adsorbatc 
induced hole resonance with density of states, po- The higher 
barrier for tunneling into the hole resonance compared to tun- 
neling into Fermi-level states, means that only a fraction of 
the total tunnel current will pass through the hole resonance, 
(b) Schematic illustration of relative energy dependent proba- 
bilities, pn(e), for inelastic hole tunneling with energy transfer 
nhtoo to the adsorbate. 



FIG. 2. (a) A single line of desorbed H from the 
Si(100)-H(2xl) surface as a result of a line scan at —7 V 
and 3.0 nA. The picture is obtained at —2 V and 0.2 nA, and 
bright regions are due to increased density of states of H-free 
Si atoms, (b) The experimental corrugation(solid line) along 
the line in (a) compared with the theoretical corrugation of a 
missing H defect including electric field effects (dashed line) 
and without electric field effects (dotted line). Results in the 
range — 7A < x < 9A are calculated using a c(4x4) cell with 
one missing H atom and outside this range using a (2x 1) cell. 
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FIG. 3. Current, Ides, as a function of sample bias, Vb, 
for constant desorption rate R = 4 s _1 . Circles show experi- 
mental results, and the solid line shows the theoretical result. 
The dashed line shows the current as a function of sample 
bias for constant electric field E = 0.8 V/A in the tunnel gap. 
Inset shows R(I) for Vb = —1 V(triangles) and —5 V(crosses), 
and from least-squares fits of R oc I a to the data we obtain 
a = 5.7 ± 0.7 (-7 V) and 6.3 ± 1.3 (-5 V). Lines show the- 
oretical calculations and have exponents a = 6.4 (—7 V) and 
7.1 (-5 V). 
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FIG. 4. a) The par- 

tial DOS, n^aie) — \(5a\fi)\' 2 5(e — e M ), of the 5a wave 
function of a Si-H molecule. The dotted line shows the pro- 
jection onto all Si-H molecular wave functions. Crosses show 
values of n 5CT (e M )/(5o-|^)| 2 . 
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